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ABSTRACT 



Aims. Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler's 
SNR and, in particular, to predict the y-ray spectrum expected from this SNR. 

Methods. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) 
explosion energy Es„ are used to constrain the astronomical and particle acceleration parameters of the system. 

Results. Under the assumption that Kepler's SN is a type la SN we determine for any given explosion energy Ef,„ and source distance d the 
mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to 
make predictions for the expected y-ray flux. Exploring the expected distance range we find that for a typical explosion energy = 10''' erg 
the expected energy flux of TeV y-rays varies from 2 x 10"" to 10"'^ erg/(cm^s) when the distance changes from d = 3.4 kpc to 7 kpc. In all 
cases the y-ray emission is dominated by /r^-decay y-rays due to nuclear CRs. Therefore Kepler's SNR represents a very promising target for 
instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of y-rays would mean that the actual source distance is larger than 7 kpc. 

Key words, cosmic rays - acceleration of particles ~ supernovae: general - ISM: individual objects: Kepler's SNR ~ radiation mechanisms: 
non-thermal - gamma-rays: theory 



1. Introduction 

Kepler's supernova remnant (SNR) (G4.5H-6.8) is the result of 
a bright supernova (SN) in our Galaxy that exploded in 1604. 
This SNR has been extensively observed throughout the elec- 
tromagnetic spectrum (for a recent review, see Blair 2005, and 
references therein). At the same time the type of Kepler's SN 
has been debated over the years. Initially it was considered a 
type la SN, based on a study of the historical light curve of 
the SN (Baade 1943). In addition its location well above the 
Galactic plane would be unexpected for a massive progeni- 
tor star. More recently it was a rgued that the Ught curv e does 
not contradict a type II-L SN (Dog gett & Bran ch 1985), and 
lBandieral ( fl987l) proposed a bow-shock model in which a mas- 
sive star, ejected from the Galactic plane, exploded into its own 
circumstellar medium. According to this picture the remnant 
now interacts with the dense bow-shock shell that was pro- 
duced in the interaction of the proge nitor's stellar wind with 
the interstellar medium. Subsequently Bo rkowski et al ] (ll992h 
calculated a detailed model of the SNR dynamics. It was ar- 
gued that the observed morphology is in good agreement with 
the bow-shock model. However, the thermal X-ray spectra, ob- 
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tained m ore recently with AS CA jKinugas a & Tsunemi|l99S|). 
Chandra (' Hwang et a l."2000') and XMM ('Cassam-Chenai et al.l 
■ 2004 ) ■ and corresponding theoretical modeling ( Badenes et alJ 
I2OO5I) . favor a type la event. We take this as our starting point. 

The similarity of their light curves and spectra favor the 
idea that SNe of type la form a rather homogeneous class with 
small variations in absolute brightness. This leads to the as- 
sumption that the corresponding mechanical explosion energy 
E^n has also a very narrow distribution. Despite the consider- 
able progress in modeling of type la SN explosions, up to now 
there is no consensus about the mean value of the explosion en- 
ergy. Within the so-called delayed-detonation model a typical 
range E^n - (1.3 - 1.6) x 10^' erg was obtained (Gamezo et al.| 
12005). This agrees w ith the range obtained from several objects 
llWheeleretalJl99i. 

The deflagration model has resulted in consider ably lower 
mean energy releases Esn = (0.4 - 0.6) x 10^' erg jKhokhlo'vl 
|200(A iReinecke et alJEool . The latter calculations however 
appear to underesti mate to s ome extent the value of the explo- 
sion energy ( Reinec ke et aI.l Eo02: Hillebrandt 2004), which is 
now expected to be rather in the range - (0.8 - 1) x 10^' erg 
in this approach (Hillebrandt, private communication). In this 
situation we use below the value E^a = 10^' erg as a typical 
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explosion energy for type la events. Since the value of 
strongly influences the SNR dynamics and in particular the 
expected y-ray flux, we explore the range £5,1 - (0.5 - 2) x 
10^' erg, in order to demonstrate the sensitivity of the final re- 
sults to the value of E^n- 

The most recent radio study of the distance to the SNR by 
[Revnoso & Goss ( 1999) leads to a lower limit of 4.8 + 1.4 kpc 
and an upper limit of 6.4 kpc. Therefore we explore be- 
low the range d - 3.4-7 kpc. For any given pair of val- 
ues Esn and d we find the density of the ambient interstellar 
medium (ISM) from a fit to the o bserved SNR size and ex- 
pansion speed jPickel et al.lll988\ This makes it possible to 
make quite definite predictions for the cosmic ray (CR) and 
■y-ray production in this SNR. We apply here the nonlinear ki- 
netic theory of CR acceleration in SNRs teerezhko et al.'l996^, 
Iperezhko & Volk 1997). As was successfully done for the rem- 
nants SN 1006, Cas A and Tycho 's SNR ( Berezhko et al. 20021 
l2003albt IVolk et al.ll2002i l2005t iBerezhko & Volk..2004l) . we 
use observations of the nonthermal radio and X-ray emission 
spectra to constrain the astronomical parameters as well as the 
particle acceleration parameters of the system, such as the in- 
terior magnetic field strength and the CR injection rates. We 
show that in all the cases considered the expected y-ray flux 
is at a detectable level if the source distance is not larger than 
7 kpc. Therefore the detection of TeV y-rays from Kepler's 
SNR (that is expected in the near future) will enable us to de- 
termine the SN explosion energy and the source distance. 

2. Model 

Hydrodynamically, a SN explosion ejects a shell of matter with 
total energy E^n and mass Mej. During an initial period the 
shell material has a broad distribution in velocity v. The fastest 
part of this ejecta distribution is describe d by a power law 
dMf .Jdv PC v^~ ^ , follo wing the treatments bv ljones et alJ(ll98ll) 
and IChevalieJ (Il982h . The interaction of the ejecta with the 
ISM - assumed to be uniform - creates a strong shock which 
heats the thermal gas and accelerates particles diffusively to a 
nonthermal CR component of comparable energy density. For 
SN la explosions we use the values Mgj - 1.4Mq and k - 1 . 
Note that since Kepler's SNR is already in the Sedov evolu- 
tionary phase its dynamics are not very sensitive to the ejecta 
parameters Mgj and k. 

Our nonlinear kinetic theory for this process is based on 
a fully time-dependent, spherically symmetric solution of the 
CR transport equations, coupled nonlinearly to the gas dynamic 
equations for the thermal gas component. Since all relevant 
equations, initial and boundary conditions for this model have 
already been described in detail in the above papers, we do not 
present them here and only briefly discuss the most important 
aspects below. 

The coupling between the ionized thermal gas (plasma) and 
the energetic particles occurs primarily through magnetic field 
fluctuations carried by the plasma which scatter energetic par- 
ticles in pitch angle. In shock waves this allows the diffusive 
acceleration of particles. The plasma physics of the field fluctu- 
ations is not yet worked out in full detail. Nevertheless, the ac- 
celerating CRs strongly excite magnetic fluctuations upstream 



of the outer SN shock in the form of Alfven waves (e.g. iBel^ 
ll978HBlandford & Ostrike3ll978l) as a result of the so-called 
streaming instability. In quasilinear approximation the wave 
amplitudes 6B grow to very high values 6B > B, w here B 
is the average field strength dMcKenzie & Volklll982h . Since 
such fluctuations scatter CRs extremely effectively, the CR 
diffusion coefficient can therefore be assumed to be as small 
as its asymptotic limit K{p) - K(mc)(p/mc), the Bohm limit, 
where K(mc) = mc^/(3eB), e and m are the particle charge and 
mass, p denotes the particle momentum, B is the magnetic field 
strength, and c is the speed of light. 

If BisM is the pre-existing field in the surrounding medium, 
then the strong streaming instability would suggest that the 
instability growth is restricted by nonlinear mechanisms to a 
level 6B ~ Bism, any further turbulent energy being dissi- 
pated into the thermal gas ( Volk & McKenzie 1981). Earlier 
attempts to give a full nonlinear description o f the magnetic 
field evolution through a numerical simulation JLucek & Bel^ 
Bell & Lucek 2001) concluded that a considerable am- 
plification of the "average" magnetic field in the smooth shock 
precursor - produced by the finite CR pressure - should occur. 
It was expected that a non-negligible fraction of the shock ram 
pressure poV^ (po is the ISM de nsity) is con verted into mag- 
netic field energy. Subsequently. i Belll il2004h argued that this 
amplification is the result of a nonresonant instability, giving 
rise to the effective upstream field Bq > Bism, on top of which 
the Alfven waves grow to amplitudes 6B ~ Bq. The Bohm limit 
and the gas heating are then to be calculated with B = Bq. 

In our analyses of the synchrotron spectrum of SN 1006, 
Tycho's SNR, and Cas A, such magnetic field amplifications 
were indeed found; they can only be produced as a nonlinear 
effect by a very effectively accelerated nuclear CR component. 
Its energy density, consistent with all existing data, is so high 
that it is able to strongly excite magnetohydrodynamic fluctu- 
ations, and thus to amplify the upstream magnetic field Bism 
to an effective field Bq > Bism, at the same time permitting 
efficient CR scattering on all scales, reaching the Bohm limit. 
The same large effective magnetic field turns out to be required, 
within the errors, by the comparison of this self-consistent the- 
ory with the field strength derived from the morphology of the 
observed X-ray synchrotron emission, in particular its spatial 
fine structure. We therefore also allow for the possibility of an 
amplified field and assume the corresponding Bohm limit for 
the particle diffusion coefficient. 

We describe the number of suprathermal protons injected 
into the acceleration process by a dimensionless injection pa- 
rameter ?/ «; 1 which is a fixed fraction of the gas particles 
entering the shock front. For simplicity we assume that the in- 
jected particles have a v elocity four times higher than the post- 
shock sound speed (see iMalkov & Dr'i^bOOlL for a review). 
It is expected that ion injection is quite efficient at the quasi- 
parallel portion of the shock surf ace, where it is c haracterized 
by values 77 = 10""* to 10"^^ (see IVolk et al.ll2003i . for details). 
Since this injection is expected to be strongly suppressed at the 
quasiperpendicular parts of the shock, one should renormaUze 
the results for the nucleonic spectrum which were calculated 
within the spherically symmetric model. The deviation from 
spherical symmetry, introduced by the divergence-free nature 
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of the magnetic field in the actual SNR, can be approximately 
taken into account by a renormalization factor < 1 , roughly 
/re = 0.15 to 0.25, which diminishes the nucleonic CR produc- 
tion efficiency calculated in the spherical model, and all effects 
associated with it. 

We assume also that electrons are injected into the acceler- 
ation process at the shock front. Formally their injection mo- 
mentum is taken to be the same as that of the protons. Since 
the physical processes of electron injection are still poorly 
known, we choose the electron injection rate such that the elec- 
tron:proton ratio K^-p (which we define as the ratio of their dis- 
tribution functions fe(p)/f(p) at all rigidities where the pro- 
tons are already relativistic and the electrons have not been yet 
cooled radiatively) is a constant to be determined from the syn- 
chrotron observations. 

The electron distribution function fe;(p) deviates only at 
sufficiently large momenta from the relation fe(p) = K^^fip) 
as a result of synchrotron losses, which are taken into account 
by supplementing the ordinary diffusive transport equation by 
a radiative loss term. 

From the point of view of injection/acceleration theory, we 
must treat K^^, together with Bo and 77, as a theoretically until 
now not very precisely calculable parameter, to be quantita- 
tively determined by comparison with the synchrotron obser- 
vations (see Volk 2003; Berezhko 2 0Q5t for reviews) 

The solution of the dynamic equations at each instant of 
time yields the CR spectrum and the spatial distributions of 
CRs and thermal gas. This allows the calculation of the ex- 
pected fluxes of nonthermal emission produced by the acceler- 
ated CRs. 

3. Results and discussion 

As a starting case we consider the typical value Esa = 10^' erg 
for the explosion energy and the conventional distance d - 
4.8 kpc. The results of our calculations are presented in Figs.^ 
|2] The hydrogen number density A^h - 3 cm"-^, which deter- 
mines the ISM density pq - \ Am^Nn, was chosen to fit the size 
/?s and the expansion speed Vj at the present age - 400 yr. 
Such a value of A^h is clearly a spatial average over the actual 
density structure. Optical Hor + [N II] data show knots and fil- 
aments, als o influenced by variable foreground extinction (e.g. 
Blair 2005, for a review). For the acceleration of the very high 
energy particles which give rise to the hard X-ray synchrotron 
and TeV y-ray emissions such small-scale inhomogeneities, 
however, play essentially no role. 

As shown in Fig. ^ the current evolutionary phase of 
Kepler's SNR corresponds to the Sedov phase. The adopted 
proton injection rate rj - 1 .5 x 10"^ leads to a significant shock 
modification, characterized by a total shock compression ra- 
tio <j ~ 6.9 and a subshock compression ratio as ~ 2.9 (see 
Fig. [fl)). Such a shock modification is needed to fit the ob- 
served steep radio spectrum and the smooth connection with 
its X-ray part (see below). 

About 10% of the explosion energy has been transfered into 
CR energy up to now, which means that the CR energy content 
is - O.lE^n- 
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Fig. 1. (a) Shock radius R^, contact discontinuity radius R^, 
shock speed V^, and contact discontinuity speed Vc, as func- 
tions of time since explosion; (b) total shock (cr) and subshock 
(cTs) compression ratios. The dotted vertical line marks the cur- 
rent epoch tc- The observed mean size and speed of the sh ock, 
as determined by radio measurements jPickel et al.lll988i) . are 
shown as well. 

The calculated synchrotron flux is shown in Fig. |2] to- 
gether with the observed values at radio and X-ray frequen- 
cies. At radio frequencies the synchrot ron spec trum 5v oc v'" 
has spectral index a - 0.71 (DeLane v et alj E002). It devi- 
ates significantly from the value a = 0.5 that corresponds to 
an unmodified strong shock. The adopted proton injection rate 
T] - 1 .5 X 10^^ gives the required shock modification. The elec- 
tron-to-proton ratio Kgp = 1.3 X 10"^ and an interior magnetic 
field strength - 480 fxG give a good fit for the experi- 
mental data in the radio and X-ray ranges. Note that the in- 
terior magnetic field = 480 fiG, derived here from the fit 
of the overall synchrotron spectrum, is higher than the value 
B'^ = 215 yuG, determined from the observed spatia l fine struc- 
ture of the synchrotron emission (Volk et al. 2005). However, 
if we use the thinne st X-ray radial profile observed by Chandra 
(Bamba et alJl2005h . which has an angular width Aif/ - 2.1", 
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Fig. 2. Calculated energy flux of synchrotron emission as a 
function of frequency for the same ca s e as in Fig. [2 The ob- 
served non-thermal X-ra vlAllen et all ( 1 19991) and radio emis- 
sion |^ynold^^^lliso3Jl993flu^^^ are also shown. 



we obtain B'^ — 340 yuG in rough agreement with the value de- 
rived from the fit of the overall synchrotron spectrum. Such a 
high interior magnetic field is the result of field amplification 
bv the nonlinear CR backreaction on the acceleration process 
llBell & Lucekll200l[ iBelf 2004). It was recently estabhshed 
that such strong field amplification takes place in all young 
Galactic SNRs which have known filanientary structures in the 
nonthermal X-ray emission dVolk et alJ2005l) . 

In Fig. 13 we present the gamma-ray spectrum of Kepler's 
SNR, expected at the current epoch. It is mainly produced by 
the CR proton component in hadronic collisions with back- 
ground gas nuclei, leading to 7r°-production and subsequent 
decay into two gamma-quanta. This so-called hadronic y-ray 
component exceeds the leptonic y-ray component due to the 
Inverse Compton (IC) scattering off the cosmic microwave 
background by more than a factor of 10^. The integral gamma- 
ray spectrum is expected to be very hard, Fy oc e^*^^, within the 
energy range from 1 GeV to almost 10 TeV. At ey = 1 TeV 
eyFy ^ 5 X 10"'^ ergAcm^s) for E^^ = lO'*' erg. Since 
the SN explosion energy is not exactly known, we present in 
Fig. 13 also the results calculated for the three other values 
£^sn/(10^' erg)=0.5, 1.5 and 2. We note that even at the low- 
est explosion energy - 0.5 x 10^' erg considered here, 
the expected y-ray flux exceeds the sensitivity of the GLAST 
instrument at GeV energies and of the H.E.S.S. instrument 
at TeV energies. At TeV-energies the expected energy flux is 
eyFy « 10"'^ erg/(cm^s) in the case /Ssn - 0.5 x 10^' erg and 
an order of magnitude higher for E^^ = 2x10^' erg. 

Qualitatively, the dependence of the expected y-ray flux on 
the explosion energy can be understood if one takes into ac- 
count that the three relevant parameters - distance d, explosion 
energy and ISM density A^h - are connected by the relation 
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Fig. 3. Total (7r"-decay + IC) integral y-ray energy fluxes as 
a function of y-ray energy for the source distance d - 4.8 kpc 
and four values of the SN explosion energy iisn/ClO^' erg)=0.5, 
1, 1.5, 2. For comparison, the respective sensitivities for a 5cr 
detection in one year with GLAST (Weekes 2003), and in 50 
hours at 20° zenith angl e for a ph oton index 2.6 (as for the Crab 
Nebula) with H.E.S.S. (lFunkl2005i) . are shown. 



because in the Sedov phase the SNR radius cc d v& deter- 
mined by the relation oc {E^JN-aS^^^t^^'^ . Taking into account 
the dependence of the Tr'^-decay y-ray flux on the relevant pa- 
rameters 



(2) 



together with the fact that in the Sedov phase the ratio E^^jE^^ 
is a weakly dependent function of time, we have 



Fy CC Eljd\ 



(3) 



(1) 



For a given source distance d this gives Fy oc E^^, in rough 
agreement with the numerical results presented in Fig.|3 

Note that magnetic field strengths vary very little with the 
explosion energy: for the interval Esn - (0.5 - 2) x 10^' erg 
considered, we have = 440-500 yuG and B'^ = 346-337 //G, 
whereas the electron:proton ratio /Tep - (2.8 -0.7) x 10""* varies 
roughly inversely proportional to the energy E^„. 

Since the source distance is not known very well, we per- 
formed our calculations for a range of distances li = 3.4-7 kpc 
in a similar way as it was done above for d - 4.8 kpc. In each 
case we achieve the same quality of fit of the observed SNR 
size, its expansion speed and the overall synchrotron emission 
spectrum. Therefore we present in Fig. |4]only the results of 
the y-ray energy fluxes expected for the SN explosion energy 
/isn = lO''' erg and for four different distances from the range 
d - 3.4-7 kpc. It can be seen from Fig.|4]that Kepler's SNR is 
expected to be as bright a TeV y-ray source as the Crab Nebula 
if the distance is as small as = 3.4 kpc. The expected y-ray 
flux goes down with increasing distance and comes to the min- 
imum observable H.E.S.S. flux if the distance becomes as large 
as 7 kpc. 
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Fig. 4. Same as in Fig.|3j but different curves correspond now 
to different source distances d - 3.4, 4.8, 6.4, 7 kpc for the 
explosion energy 10^' erg. 

According to Eq. 0, for a fixed value of the SN explosion 
energy /Ssn the expected -y-ray flux is 

Fy CC (4) 

in good agreement with the results presented in Fig.0] 

We note he re that according to the so-called delayed det- 
onation model dOamezo et al.ll2005l) the value E^n — 10^' erg 
is even below the lower end of the expected range of type la 
SN explosion energies. From this point of view the results pre- 
sented in Fig. 0] most probably give a lower limit for the ex- 
pected y-ray flux for any given distance d, and therefore these 
results convince us that Kepler's SNR is a potential high-energy 
y-ray source. 

We emphasize that an increase of the assumed source dis- 
tance leads to an increase of the linear width L oc dAifr of the 
filaments observed in nonthermal X-rays. This leads in turn to a 
decrease of the interior magnetic field strength B'^ oc L"^^^, de- 
termined from the filament size. On the other hand, the linear 
shock size R^^ cc d and the shock speed Vj oc become larger 
for larger d and this makes the CR acceleration process faster. 
Therefore a higher magnetic field with a correspondingly 
more intense synchrotron cooling effect of CR electrons is re- 
quired in order to provide a smooth cutoff of the synchrotron 
spectrum 5y(v) consistent with the existing nonthermal X-ray 
measurements (see Fig.|3. As a result, at li = 7 kpc we have 
Bd = 534 fiG, whereas B'^ - 295 fiG is already considerably 
lower This means that a source distance d x 1 kpc or larger 
has to be considered as inappropriate if the radial X-ray profile 
width Al// - 2.1" observed by Chandra (Bamba et al...2005i) is 
indeed the smallest one in this source. Then Kepler's SNR has 
to be considered as a promising potential source of high-energy 
y-rays . We note also that the electron:proton ratio Agp is quite 
insensitive to the distance d. 

It is clear that the asymmetry of the X-ray emission is ev- 
idence that the surrounding ISM is inhomogeneous on a scale 
comparable with the SNR size (e.g. lBandiera 19871 Our spher- 



ically symmetrical approach is not able to take this into ac- 
count. However, since we used the SNR size Rs and speed Vs 
as averages over the remnant, we believe that the ISM density 
adopted for each pair of values d and E^n also represents the 
ambient density averaged over the remnant and thus furnishes 
a reliable estimate for the expected total gamma-ray emission. 

4. Summary 

Our consideration of CR acceleration and nonthermal emis- 
sion in Kepler's SNR demonstrates that a spherically symmet- 
ric nonlinear kinetic theory reproduces the SNR dynamics and 
the properties of its nonthermal radiation in a very satisfactory 
way when Kepler's SN is treated as a type la SN. 

The y-ray energy flux expected at TeV energies is eyFy « 
(3 - 5) X 10"'^ erg/(cm^s) if the distance is as small as d - 
4.8 kpc. The flux is expected to be in a detectable range eyFy > 
10"'^' erg/(cm^s) at TeV energies if the distance does not ex- 
ceed 7 kpc. If the upper limit for th e source distance is indeed 
d - 6.4 kpc jRevnoso & Gossll999 ) - a conclusion that is con- 
firmed by the consistency check of the interior magnetic field 
values obtained by two independent methods - we conclude 
that Kepler's SNR is a potentially bright y-ray source in the 
sky. 

Acknowledgements. This work has been supported in part by the 
Russian Foundation for Basic Research (grant 03-02-16524) and LSS 
422.2003.2. EGB and LTK acknowledge the hospitality of the Max- 
Planck-Institut fiir Kernphysik, where part of this work was carried 
out. The authors thank W. Hillebrandt and G. Puhlhofer for discus- 
sions on type la SN explosions and the specific properties of Kepler's 
SN. 

References 

Allen, G. E., Gotthelf, E. V., & Petre, R. 1999, Proc. 26th 

ICRC, Salt Lake City, 3, 480 
Baade, W. 1943, ApJ, 97, 119 

Badenes, C, Borkowski, K. J., & Bravo, E. 2005, ApJ, 624, 
198 

Bamba, A., Yamazaki, R., Yoshida, T., Terasawa, T., & 

Koyama, K. 2005, ApJ, 621, 793 
Bandiera, R. 1987, ApJ, 319, 885 
Bell, A. R. 1978, MNRAS, 182, 147 
Bell, A. R. 2004, MNRAS, 353, 550 
Bell, A. R. & Lucek, S. G. 2001, MNRAS, 321, 433 
Berezhko, E. G. 2005, Advances in Space Research, 35, 1031 
Berezhko, E. G., Elshin, V. K., & Ksenofontov, L. T 1996, J. 

Exp. Theor. Phys., 82, 1 
Berezhko, E. G., Ksenofontov, L. T, & Volk, H. J. 2002, A&A, 

395, 943 

Berezhko, E. G., Ksenofontov, L. T., & Volk, H. J. 2003a, 
A&A, 412, LI 1 

Berezhko, E. G., Puhlhofer, G., & Volk, H. J. 2003b, A&A, 
400, 971 

Berezhko, E. G. & Volk, H. J. 1997, Astroparticle Physics, 7, 
183 

Berezhko, E. G. & Volk, H. J. 2004, Astron.Astrophys., 419, 
L27 



6 



E.G.Berezhko et al.: Gamma-ray emission expected from Kepler's SNR 



Blair, W. P. 2005, in ASP Conf. Ser. 342: 1604-2004: 
Supernovae as Cosmological Lighthouses, 416 (also 
astro-ph/0410081l 
Blandford, R. D. & Ostriker, J. P 1978, ApJ, 221, L29 
Borkowski, K. J., Blondin, J. M., & Sarazin, C. L. 1992, ApJ, 
400, 222 

Cassam-Chenai, G., Decourchelle, A., Ballet, J., et al. 2004, 

A&A, 414, 545 
Chevalier, R. A. 1982, ApJ, 258, 790 

DeLaney, T., Koralesky, B., Rudnick, L., & Dickel, J. R. 2002, 
ApJ, 580,914 

Dickel, J. R., Sault, R., Arendt, R. G., Korista, K. T., & Matsui, 

Y. 1988, ApJ, 330, 254 
Doggett, J. B. & Branch, D. 1985, AJ, 90, 2303 
Funk, S. 2005, in PhD - thesis, Univ. of Heidelberg, Germany, 

71 

Gamezo, V. N., Khokhlov, A. M., & Oran, E. S. 2005, ApJ, 
623, 337 

Hillebrandt, W. 2004, New A Rev., 48, 615 

Hwang, U., Holt, S. S., Peti-e, R., Szymkowiak, A. E., & 

Borkowski, K. J. 2000, BAAS, 32, 1236 
Jones, E. M., Smith, B. W., & Straka, W. C. 1981, ApJ, 249, 

185 

Khokhlov, A. M. 2000, arXiv: astro-ph/0008463J 
Kinugasa, K. & Tsunemi, H. 1999, PAS J, 51, 239 
Lucek, S. G. & Bell, A. R. 2000, MNRAS, 314, 65 
Malkov, M. A. & Drury, L. O. 2001, Reports of Progress in 

Physics, 64, 429 
McKenzie, J. F. & Volk, H. J. 1982, A&A, 116, 191 
Reinecke, M., Hillebrandt, W., & Niemeyer, J. C. 2002, A&A, 

386, 936 

Reynolds, S. P & Ellison, D. C. 1992, ApJ, 399, L75 
Reynoso, E. M. & Goss, W. M. 1999, AJ, 118, 926 
Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2003, A&A, 
409, 563 

Volk, J. H. & McKenzie, F J. 1981, in Proc. 17th ICRC, Paris, 
246-249 

Volk, H. J. 2003, in Proc. 28th ICRC (Tsukuba), ed. T. Kajita, 
Y. Asaoka, A. Kawachi, Y. Matsubara, & M. Sasaki (Tokyo, 
Japan: Universal Academy Press Inc.), Vol. 8, 29^8 (also 
astro-ph/0312585l 

Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2005, A&A, 
433, 229 

Volk, H. J., Berezhko, E. G., Ksenofontov, L. T., & Rowell, 

G. P 2002, A&A, 396, 649 
Weekes, T. C. 2003, in Proc. 28th ICRC (Tsukuba), ed. T 

Kajita, Y. Asaoka, A. Kawachi, Y. Matsubara, & M. Sasaki 

(Tokyo, Japan: Universal Academy Press Inc.), Vol. 8, 3-15 
Wheeler, J. C, Harkness, R. P, Khokhlov, A. M., & Hoflich, 

P 1995, Phys. Rep., 256, 211 



